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Abstract. We obtain a formula for the 71-dimensional distributions of the Airyi process 
. in terms of a Fredholm determinant on I/^(R), as opposed to the standard formula which 

involves extended kernels, on L^{{1, . . . ,n} x R). The formula is analogous to an earlier 
^— ^ I formula of Prahofer and Spohn [PS02] for the Airy2 process. Using this formula we are 

■ able to prove that the Airyi process is Holder continuous with exponent i— and that it 

^ ■ fluctuates locally like a Brownian motion. We also explain how the same methods can be 

used to obtain the analogous results for the Airy2 process. As a consequence of these two 
results, we derive a formula for the continuum statistics of the Airyi process, analogous 



< 



' to that obtained in [CQRll] for the Airy2 process. 



! 1. Introduction and Main Results 

CIh; 

1.1. General background. The Airy processes are stochastic processes which are ex- 
pected to govern the asymptotic spatial fluctuations in a wide variety of random growth 
models on a one dimensional substrate, top lines of non-intersecting random walks and free 
energies of directed random polymers in 1 -|- 1 dimensions (all belonging to the Kardar- 
Parisi-Zhang, or KPZ, universality class [KPZ86]). They are non-Markovian and are de- 
, fined in terms of their finite-dimensional distributions, which are given by determinantal 

Q>^ ■ formulas. These formulas, which have been derived by asymptotic analysis of exact for- 

■ mulas in special discrete models such as the totally asymmetric simple exclusion process 

I and the polynuclear growth model, give the n-dimensional distributions in terms of Fred- 

holm determinants of extended kernels, on L^({1, . . . , n} x M). The exact results are then 
^ . conjecturally extrapolated to more general processes in the universality class which do not 

^ I possess the same exact solvability. 

The particular Airy process arising in each case depends on the initial data, and this 
picks out a number of KPZ sub-universality classes. For reasons of scaling invariance, there 
are three special pure initial data classes: narrow wedge, flat, and equilibrium. Narrow 
^ ■ wedge corresponds to point-to-point polymers, or growth models where the exponential of 

5^ . the height is initially a Dirac delta. Physically, one starts with curved, or droplet, initial 

data. After some time t, the height looks like a parabola in space, corresponding to the 
deterministic evolution, on top of which is approximately an Airy2 process [PS02] with 
amplitude t^^^ and varying on a spatial scale of t^^^. Flat corresponds to point-to- line 
polymers, or growth models with constant initial data. At time t, one sees spatially the 
Airyi process [Sas05], again with size t^^^ and varying on spatial scale t'^^^. Equilibrium 
corresponds to growth models starting from equilibrium, which in the KPZ universality 
class means approximately a two-sided Brownian motion. At a later time one sees spatially 
the Airystat process [BFPIO]. Note that all these descriptions are modulo a global height 
shift which is non-trivial itself, and can be very large compared to the scales on which these 
fluctuations are observed. 

There are also three other basic mixed initial data, corresponding to starting with one 
of the basic three geometries to the left of the origin and another one to the right. The 
resulting spatial fluctations are still of size t^^^ and on a spatial scale of t^/^, with non- 
homogeneous crossover Airy processes Airy2_>i [BFS08], Airyi_j>stat [BFS09] and Airy2_s>stat 
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[SI04; CFPIO], the names being self-explanatory. Of course, there will be other less com- 
monly seen sub-universality classes, but these six are the basic ones, and, interestingly, all 
have determinantal finite-dimensional distributions. 

Although the determinantal formulas arise naturally in deriving the finite-dimensional 
distributions from the special solvable discrete models, they are cumbersome for the anal- 
ysis of properties of these processes involving short range scales. For example, one would 
expect to be able prove the pathwise continuity directly by just checking the Kolmogorov 
continuity criterion using the determinantal formula for the two point distributions with 
extended kernel on L^({1,2} x M). This turned out to be surprisingly difficult, and has 
been an open problem since the processes were introduced. For the Airy2 process, which is 
in some sense the most basic one, what was done historically was to study the probability 
measure on the point processes obtained by sampling the Airy line ensemble at a finite 
set of times. Prahofer and Spohn [PS02] proved the continuity of the Airy line ensemble 
as a point process, from which the continuity of the top line, the Airy2 process, would 
follow if one knew that the points came from a non-intersecting line ensemble. However, 
this was not known at the time (though it is now, see [CHll]). Johansson [Joh03] proved 
the tightness of an approximating line ensemble (the multilayer PNG model), which in 
particular implied the continuity of the Airy2 process. 

On the other hand, the other processes do not arise easily as top lines of line ensembles. 
For example, for the Airyi process, which will be our main example in this article, even 
the continuity remained open. 

One also hopes to study variational problems involving the Airy processes. These arise 
naturally. A well-known example is the famous result of Johansson [Joh03] that the 
supremum of the Airy2 process minus a parabola has the Tracy- Widom GOE distribu- 
tion [TW96]. There is also a generalization of this [QRllb] that the same supremum on a 
half-line is given by the one point marginal of the Airy2_).i process. Variational problems 
naturally involve infinitely many spatial points, so formulas giving the distribution of n 
sample points in terms of determinants of extended kernels on L^({1, . . . , n} x M) are not 
a good tool. In [CQRll] we introduced a continuum formula for the Airy2 process, which 
gives the probability that the process lies below a given function on an arbitrary finite in- 
terval, in terms of a Fredholm determinant of the solution operator of a certain boundary 
value problem. The formula is obtained as a fine mesh limit of an older formula of Prahofer 
and Spohn [PS02] for the n-dimensional distributions (see (1.6) below). The advantage of 
the alternative formula for variational analysis is that its complexity is no longer diverging 
with the number of spatial points. Using this formula, we were able to give a direct proof 
of Johansson's result [Joh03], study the half line version [QRllb], and derive an exact 
formula for the probability density of the argmax of the Airy2 process minus a parabola, 
the polymer endpoint distribution [MQRll]. 

In this article we will obtain analogous discrete and continuum formulas for the Airyi 
process, and use them to prove directly that it is Holder ^ — 6 continuous for any 6 > 0. 
This regularity of Airyi is expected from the fact that the process is believed to look locally 
like a Brownian motion. In fact, we will show in this direction, using the alternative deter- 
minantal formula, that the finite dimensional distributions of the Airyi process converge 
under diffusive scaling to those of a Brownian motion. 

Note that the existence of formulas for the Airyi process involving boundary value 
operators is to some extent surprising. In the case of the Airy2 process, which is the 
limit of the rescaled top line in a system of non-intersecting Brownian motions (Dyson's 
Brownian motion for the Gaussian Unitary Ensemble), the formula can be seen as a certain 
extension of the Karlin-McGregor formula (see [BCR12]). On the other hand, there is no 
known analogous construction of the Airyi process (see in particular [BFP08]), for which 
the associated determinantal process is signed (see [BFP07]), and thus it is not at all 
apparent where formulas like (1.7) or (1.14) below are coming from. 
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1.2. Statement of the results. Now we turn to a precise description of the Airyi process, 
which win be our main object of study. It was first derived by Sasamoto [Sas05] (see also 
[BFPS07; BFP07]) by asymptotic analysis of exact formulas for TASEP with periodic 
initial data. It is a stationary process defined through its finite-dimensional distributions, 
given by a determinantal formula: for xi, . . . ,Xn G M and ti < • • • < t„ in M, 

(1.1) ¥{Ai{h)<xi,...,Ai{tn)<Xn) = det{I-f/^KrH'/^)L2^{t,,..^^^^^^ 

where we have counting measure on {ti, . . . and Lebesgue measure on M, f is defined 

on X M by f{tj,x) = la;e(xj,oo) and 

(1.2) irr*(t,x;t',x') = , ^ expf- ^f'~'^n 

+ Ai(a; + x' + {t' - tf) exp (^t' - t){x + x') + ^(i' - • 

Here, and in everything that follows, the determinant means the Fredholm determinant in 
the Hilbert space indicated in the subscript. In particular from (1.2) and [FS05] one obtains 
that the one-point distribution of the Airyi process is given in terms of the Tracy- Widom 
largest eigenvalue distribution for the Gaussian Orthogonal Ensemble (GOE) [TW96]: 

P(^i(0) <m) =FGOE(2m). 

Note that it follows from (1.1) that Ai{t) has the same distribution as Ai{—t). 

The definition of the Airyi process is analogous to that of the Airy2 process, introduced 
by Prahofer and Spohn [PS02], whose n dimensional distributions are given by 

(1.3) P(^2(ti) <rEi,...,^2(tn) <x„) = det(/-fV2ii'-tfV2)^,^^^^ 

where the extended Airy kernel [PS02; FNH99; Mac94] Kf"^ is defined by 

, , , //o°°(iAe-^(*-*')Ai(x + A)Ai(x' + A), if t > t' 
^ ^' ' '^"\-/\dAe-^(*-*')Ai(x + A)Ai(x' + A), ift<t',- 

The analogy between the definitions becomes clearer in light of the following observa- 
tions. Letting Kx\ denote the Airy kernel 

KAi{x,y)= [ dAAi(x- A)Ai(y- A) 



and H denote the Airy Hamiltonian 

H = -A + x, 

where A = 5^ denotes the one-dimensional Laplacian, one can show (formally) that the 
extended Airy kernel can be rewritten as 

(1.4) Kl^\t,x;t',x') = -e-^''-''^" {x,x')lt>>t + e'^ KAie"''" {x,x'). 

On the other hand, as shown in Appendix A of [BFPS07], K^^^ can be expressed (formally) 
in the following alternative way: 

(1.5) i^f*(*,a;;i',x') = -e(*'-*)^(x, x')lt'>t + e-*^5oe*''^(x, x'), 
where 

Bo{x,y) =Ai{x + y). 

Note that (1.5) corresponds exactly to (1.4) after replacing H by —A and K^i by Bq. This 
particular replacement was emphasized in [Fer08]; more generally, all the extended kernels 
arising in this and related areas have an analogous structure. We stress that both (1.4) 
and (1.5) should be regarded at this point as formal identities, as it is not clear how to 
make sense of e^*^ and e*'^ for t < 0. 
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Our first result provides a new determinantal formula for the finite-dimensional distri- 
butions of the Airyi process without using extended kernels or, in other words, involving 
the Fredholm determinant of an operator acting on L^(M) instead of . . . ,t„} x M). 

For the Airy2 process such a formula was introduced by Prahofer and Spohn [PS02] as its 
original definition: 

(1.6) nA2{tl)<Xi,...,A2{tn)<Xn) 

= det(/-KAi + P.,e(*-*^)^P.,e(*-*3)H...p^^^(t„-t.)/.^ \ 

where Pa denotes projection onto the interval (—00, a]. The equivalence of (1.3) and (1.6) 
was derived in [PS02; PSll], see Remarks 2.1 and 2.2 below for a discussion about some 
technical details. Our result states that the finite-dimensional distributions of the Airyi 
process admit the same representation after replacing by — A and A'ai by Bq. 

Theorem 1. The finite- dimensional distributions of the Airyi process are given by the 
following formula: for xi, . . . , a;„ € R and ti < • • • < t„ in M, 

(1.7) F{Ai{ti)<xi,...,Ai{tn)<Xn) 

= det(/ -Bo + P,,e-(*i-*2)^/>,,e-(*^-*^)^ • • • P,„e-(*"-*^)^i?o) ^^^^^ • 

Remark 1.1. 

1. Note that, since ti < ■ ■ ■ < tn, all the heat kernels in (1.7) are well defined except 
for the first one. The same situation is present in the formula for the Airy2 process, 
as the factor e(*"~*i)^ in (1.6) is in principle ill-defined. The situation is resolved 
in that case by observing that e.i^n-ti)H -g applied after Kp^\ in (1.6), and Kp^\ is a 
projection operator onto the negative eigenspace of H. In our case the situation is 
resolved by Proposition 1.2 below. 

2. The operator 

appearing inside the determinant in (1.7) is not trace class, basically because the 
heat kernel is not even Hilbert-Schmidt. However, we will show in Proposition 
2.3 that there is a conjugate operator J = U~^JU which is trace class in L^(M), 
so the formula (1.7) should be computed as det(7 — J)i2(j^y Alternatively, this 
implies that the Fredholm determinant in (1.7) regarded as its Fredholm expansion 
series is well defined. (The same issue arises in (1.1), as Kf^^ is not trace class on 
L2({ti, ... ,tn} xR); this is resolved in Appendix A of [BFP07]). 

3. Note that the issue discussed in the last point does not arise in the formula (1.6) 
for the Airy2 process. The fact that the operator appearing in that formula is trace 
class is proved in Proposition 3.2 of [CQRll]. 

The following result shows that we are allowed to consider the operator e~*^ for t > 
as long as it is applied after Bq. 

Proposition 1.2. For fixed t,y eM let ^t,y{x) = e-2*V3-(x+y)t Ai(a; + y + t^). Then for 
all s,t > we have 

(1.8) e''^ipt,y{x) = ipt-s,yix). 

In particular, e^^(ft,y = Ai(a; + y), and as a consequence the kernel c^^^Bq is well defined 
for every t > via the formula 

(1.9) e-'^^Bo = e-2*V3-(x+,)t Ai(^ + y + t^) 

and it satisfies the group property in the sense that c^^^^^^Bq = c^^c^^Bq for all s, t € R. 
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We remark that versions of the above identities appear in earher works on the Airyi 
process, and in particular in [Sas05; BFP07; BFPS07]. Proposition 1.2 allows us to make 
sense of (1.5): since the Airyi process is stationary, by shifting ti, . . . we may assume 
that < ti < ■ ■ ■ < tn, and then all the heat kernels with a negative parameter in (1.5) 
appear applied after Bq. The same type of argument allows to make sense of (1.4), (1.6) 
and (1.7) (though see also the last paragraph of Remark 2.2). 

As we mentioned, formulas (1.6) and (1.7) are better adapted than the standard extended 
kernel formulas to short range properties of the process. As a first application we will prove 

Theorem 2. The Airyi process Ai and the Airy2 process A2 have versions with Holder 
continuous paths with exponent ^ — S for any 5 > 0. 

Recall that continuity was known for A2 but not for Ai. The Holder ^— continuity 
for A2 also follows from recent work of Corwin and Hammond [CHll]. They study the 
Airy line ensemble directly, obtaining the continuity (and Holder ^— continuity) directly 
from a certain Brownian Gibbs property. In general, all the Airy processes are supposed 
to be locally Brownian. Note that the definition of locally Brownian is not unique. For 
A2 it follows from [CHll] that it is locally absolutely continuous with respect to Brownian 
motion. Analogous results have recently become available for the solutions of the KPZ 
equation at finite times [Haill; QRlla; CH12]. For ^1 the line ensemble picture is missing 
at the present time, so a proof was lacking. As another application of the formulas, we 
prove that the Airyi process is locally Brownian in the sense that under local Brownian 
scaling, the incremental process converges to that of Brownian motion. 

Theorem 3. For any fixed s £ M, let Bs{-) be defined hy Be{t) = e~''/'^{Ai{s + et)-Ai{s)), 
t > 0. Then B^{-) converges to Brownian motion in the sense of convergence of finite 
dimensional distributions. The same holds for -Be(-) defined by Bf,{t) = Bf,{—t), t > 0. 

Note that by stationarity there is no loss of generality in taking s = in the theorem, 
while the statement about Bf,{-) follows from the statement about -Be(-) by time reversal 
invariance of Airyi. The analogue of Theorem 3 for Airy2, which follows from its local 
absolute continuity with respect to Brownian motion, was proved earlier by Hagg [Hag08], 
and can also be obtained directly by our method. We remark also that, using an analogue 
of (1.7) for the Airy2_5.i process, which will appear in upcoming work [BCR12], it should 
not be hard to adapt our proofs to show that A2^i is Holder i— continuous and is locally 
Brownian in the sense of the last result (in fact, the result of [BCR12] is more general and 
should allow one to extend our proofs to other processes). 

Going back to ^1, one can be quite precise in terms of finite dimensional distributions. 
Letting < ti < • • • < t^, we will prove that 



(1.10) ¥{Al{eti)<X + y/^yi,...,Al{£tn) <x + s/Iyn\Ai{Qi) = x) 

= lE(lB(t,)<j;„i=l,„,,„fi't,y(2;,-B(in))) hl^^{x), 



_ J^^ du e~^'-^BoiV^z + X, u) [1-B0 + A[j^^''+"-)e-^*"^i?o) (n, x) 
f^^duBo{x,u) {I - Bo + P^Bo) {u,x) 
where a[5^^^+") = P,e*i^P,,+,e(*^-*i)^ • • • e(*""*"-)^P,„+, and 

¥{Ai{0)<x,Ai{eti)<x + ^yi,...,Ai{etn)<x + ^yn) 




(1.12) 




i^GOE(22;) 



One has 



lim5fy(2;,^) = lini/if (x) = 1 
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from which it fohows from (1.10) that the finite dimensional distributions converge to 
those of Brownian motion. It would be interesting to understand the role of gly{x,z). 

Expansions gly{x, z) = 1 + e^/^fl'^^y(x, z) + 0{£) and hl y{x) = 1 + e^/^/i[^y(x) + 0{e) may- 
identify the infinitesimal increments of Ai in order to develop a stochastic calculus. 

One of course has formulas analogous to (1.10) for the Airy2 process (and, in view of 
[BCR12], other processes such as Airy2_i.i), but we do not include them here. 

Our last result, which is an application of Theorems 1 and 2, gives a determinantal 
formula for the continuum statistics of the Airyi process on a finite interval. This was 
done for the Airy2 process in [CQRll], and the same argument will allow us to take a limit 
of the formula in Theorem 1 as the size of the mesh in t goes to 0. 

Fix £ < r. Given g G H^{[£,r]) (i.e. both g and its derivative are in L'^[[£,r])), define an 
operator A^^^j acting on L^(M) as follows: A^^^j/(-) = u{r,-), where u{r,-) is the solution 
at time r of the boundary value problem 

dtu - An = for X < g{t), t £ {i, r) 

(1.13) n(Aa;) = /(x)l^<g(f) 

u{t,x) = for X > g{t). 

The fact that this problem makes sense for g G H^{[£, r]) is not hard and can be seen from 
the proof of Proposition 2.3 below (see also Proposition 3.2 of [CQRll]). 

Theorem 4. 

(1.14) ¥{Ai{t) < g{t) for t G [£,r]) = det(l - Bo + Af,^^,e-(^-^)^So 



In other words, hitting probabilities of curves by Ai can be expressed in terms of Fredholm 
determinants of the analogous hitting probabilities for Brownian motion. 

One can check easily using the Feynman-Kac formula that the kernel of A^^ has the 
following form: 

(x-y)yiir-l) 

(1.15) Af,,,,(x,,) = -^J^^=^nii)=.Mr)=y[^(') ^ on [i,r]) , 

where the probability is computed with respect to a Brownian bridge b{s) from x at time 
£ to y at time r and with diffusion coefficient 2. We remark that the kernel —Bq + 
A^^ ^^e~^'^~^^^BQ is not trace class, but as in the discrete case (see Remark 1.1) we will 
show that there is conjugate operator which is, see Proposition 2.3. 

The corresponding formula for the Airy2 process, provided in Theorem 2 of [CQRll], is 
the same as (1-14) after replacing —AhyH and Bq by K\[. The corresponding boundary 
value operator ©^^^j in that case is actually more complicated than A^^^j, as in our case 
there is no potential term in the partial differential equation in (1.13). 

Acknowledgements. Both authors were supported by the Natural Science and Engineer- 
ing Research Council of Canada, and DR was supported by a Fields-Ontario Postdoctoral 
Fellowship and by Fondecyt Grant 1120309. The authors thank Ivan Corwin and Alexei 
Borodin for interesting and useful conversations. 



2. Proof of the determinantal formula 

Throughout this section and the next we will denote by || • ||i and || • ||2 respectively the 
trace class and Hilbert-Schmidt norms of operators on L^(M) (see Section 3 of [CQRll] for 
the definitions or [Sim05] for a complete treatment). 
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Proof of Proposition 1.2. Recall that Ai(z) = (27ri) ^ Jp due 
iy, y & M} for any fixed c > 0. Then 



,U^ 



/3 uz^ where Fc = {c + 




We can compute the z integral first, which is just a Gaussian integral, to obtain 



' ^ ' 27ri r 
Shifting to n — s we get 




which proves (1.8). The remaining statements in the proposition follow directly from this 



We turn now to the proof of Theorem 1. The argument is based on the derivation of the 
equivalence of (1-3) and (1.6) for the Airy2 case given by Prolhac and Spohn [PSll], and in 
fact the algebraic procedure we will use is basically equivalent to theirs. In the case of the 
Airyi process one has to make sure throughout the proof that the algebraic manipulations 
are being done on operators which are trace class, so that the Fredholm determinants 
considered are well defined. This is done by rewriting the algebraic procedure of [PSll] so 
that in each step one can conjugate by the correct operators and check that the resulting 
conjugated operators are trace class as needed. 

Remark 2.1. Our proof of Theorem 1 can be used to complete the details and provide all 
the necessary justifications in the proof given in [PSll] for the Airy2 case. In one sense the 
argument in that case is simpler, because the kernels in (1.3) and (1.6) are already trace 
class. Nevertheless the Airy2 case presents an additional difficulty, namely that even for 
t > the operator e~*^ does not map L^(]R) into itself (note that this issue does not arise 
in the Airyi case, as e*^ is clearly a bounded operator acting on L^(M) for t > 0). We will 
explain in Remark 2.2 how this can be addressed, and in particular how the proof below 
has to be changed to provide a rigorous proof for the Airy2 case. 

Proof of Theorem 1. We will retain most of the notation of [PSll], and as in that paper 
we use sans-serif fonts (e.g. T) for operators on L^({ti, . . . x M). Such an operator 
can be regarded as an operator- valued matrix (Tij^. ^ with entries Tij € L^(M) 




acting on / G L^(M)" as (T/)j = Yl^=i'^i,jfj (o^j more precisely, as an operator acting 
on (8) L^(R)). We will use serif fonts for the matrix entries (e.g. Tjj- = T for some 
T S L^(M)). All determinants throughout this proof are computed on L^({ti, . . . , f„} x M) 
unless otherwise indicated. 

Recall from Proposition 1.2 that e*^i?o satisfies the semigroup property e^^e^^B^ = 
g(s+t)A^^ for all s, t e M. We will use this fact several times below. We will also use the 
fact that, since BQ{x,y) depends only on x + y, e*^ and Bq commute for t > 0. Finally, 
as explained after the proof of Proposition 1.2, we may (and will) assume that ti > for 
i = 1, . . . , n. 

Let K = fV2^i^ji/2^ ^j^h j^i^^ defined through (1.5) and f as in (1.1). Using the above 
interpretation K can be written as 



identity. 



□ 




(2.1) 



K = P(T^K° + T+(K° - l))P 



where 




Boh. 
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with Pa = I — Pa denoting projection onto the interval [a, oo), and T~, T"'' are lower 
triangular, respectively strictly upper triangular, and defined by 

I — e J.j>j , ^ ij — ^ ^t<3 ■ 

Observe that the all heat kernels in T"*" have positive parameters, while those in T~ have 
negative parameters but appear applied after Bq in the expression for K in (2.1), so Propo- 
sition 1.2 ensures that (2.1) makes sense. 

As we mentioned in Remark 1.1, it is proved in [BFP07] that there is an invertible 
operator V such that VKV^^ is trace class. Explicitly, V is a (diagonal) multiplication 
operator given by 

V,j=V-l,=j with Vif{x) = {l + x^)-^'f{x). 

Since VPT^PV"^ is strictly upper triangular, I + VPT'^PV"^ is invertible, and then we can 
write 

(2.2) det(l - VKV~i) = det((l + Wi)(l - (I + W^-^Ws)) 
with 

(2.3) Wi = VPT+PV-\ W2 = VP(T- +T+)K°PV-^ 
We remark that Wi is trace class by Lemma A. 2 in [BFP07]. 

Next we want to obtain an explicit expression for (I + Wi)^iW2. Observe that 

(2.4) [(I + T+)-i)]^_^. = 1,=, - e-(*-*'+^)^li=,_i, 

which can be checked directly using the semigroup property of the heat kernel. In particular 
I + is invertible, so we can write 

(2.5) (I + Wi)-^W2 = (I + Wi)"^VP(T- + T+)(l + T+)-^K°(l + T+)PV^\ 

where we have used the fact that e*^ and Bq commute for t > 0, and hence so do and 
K°. Using (2.4) we deduce that 

[(T- + T+)(l + T+)-iK°] = e-(*'-*^)^5o - e-^'^-'^-^^^e-^'^-'-'^^^ B^lj^^ 

^^■^^ = e-(*»-*^)^5ol,=i. 

Note that only the first column of this matrix has non-zero entries. 

Observe now that, since VPT+PV-^ is strictly upper triangular, we have (VPT+PV ^)"-+^ 
0, which implies that 

n 

(2.7) (I + Wi)-^ = ^(-1)*^(VPT+PV-^)^ 

A,'=0 

On the other hand by (2.6) we have for < k < n — i 



(2.8) r(VPT+PV~^)^'VP(T- + T+)(l + T+)-^K^ 



i<ai <■ ■ ■ <ak<n 

which follows from (2.6) and the definition of PT"'"P, while ior k > n — i the left side above 
equals (and the case A: = is interpreted as ViPx^e~^^^~^''-^^BQ). Replacing each factor 
Px except the first one by / — -Pr and using the semigroup property for the heat kernel we 
deduce that the last expression equals 

k 



m=0 i=bo<bi<---<bm<n 
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Summing the above expression times (—1)^ from k = to k = n — i and interchanging the 
order of summation leads to 



n—i n—i 



EE E 

m=0 k=m i=bQ<b\< - ■ -Kb^Kn 



n — I — m 
k — m 



)(-l) 



k+m 



Noting that YTJm (™ 
(2.9) 

[(I + Wi)"ivp(T- + T+)(i + j+y^v: 



and recalhng (2.7) we deduce that 

1- E 



'-m=n—i 



i=bo<bi<---<bn-i<n 



Post-multiplying by (I + T+)PV~^ we finally obtain from this and (2.5) that 
(2.10) [(I + Wi)-iW2]^ . = li=„Ki'x„e-(*"-*^)^i?oPx, V^,"' 

where we have used again the fact that e*^ commutes with Bq for t > 0. 

At this stage we can check that (I + Wi)~^W2 is trace class. In fact it is enough to check 
(see (A. 5) in [BFP07]) that each entry of this operator- valued matrix is trace class. The 
case i = n was checked in Lemma A. 3 in [BFP07], while for the case i < n we can use a 
similar strategy. Since Vi and V^^^ are multiplication operators, they commute with Pa for 
any a, and then choosing —L< minjxj, xj} we have 



[(l + Wi)-iW2] 



K;Px.. P- LRie-^'--'^^^BoP^LPx. K- 



P,T_LViRie-^'--'^^^BaVr^P^LPx, ||^ < \p^LViR^e-'^'--'^^^BQVr^P^L 

where Ri = e~(*»~*»+i)Ap^..^^e~(**+i~**+2)A • • • P^^ and we have used the first of the inequal- 
ities 



(2.11) 
with II 



l^^lll < Pllopll^lll, 



l^^lb < Pllopll^lb, 



IASIIl < Plllllslll, 



lop 



denoting the operator norm (see [Sim05]) and ||Pa 



a; 1 1 op 



1. Next we remove 

the projections P-l and think instead of the operator ViRie~^^"~^^^^ BqV~^ as acting on 
L^([— L, oo)). Using again (2.11) and the fact that the operators Vi and commute with 
Pa we have that ||Vii?iV^-'^||i is bounded by 

which is finite because each factor is so by Lemma A. 2 in [BFP07]. Since ||V^e~^*"~*^^^-Bo^^^||i 
(computed in L^([—L, oo)) is finite by Lemma A. 3 in [BFP07] we deduce by (2.11) that 

|[(l+Wi)"iW2]^J^ < \\ViR^V~%\\Vne~^'--'^^^BQVr% < ^. 

Going back to (2.2), since both Wi and (I -|- Wi)~^W2 are trace class, we have 

(2.12) det(l-VKV"i) =det(l + Wi)det(l-(l + Wi)-iW2) = det(l - (l + Wi)-iW2), 

where the second equality follows from the fact that, since Wi is strictly upper triangular, 
its only eigenvalue is 0, and thus det(l -|- Wi) = 1. Now let U be given by Ujj- = Uli=j 
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where U is the (diagonal) multiplication operator introduced in Proposition 2.3. Then x 
to (2.5) we have 

(l + Wi)-^W2 = W3W4 

with W3 = (I + Wi)-iVP(T- + T+)(l + T+)-1kOU-i and W4 = U(l + T+)PV-i. We have 
already checked that W3W4 is trace class, so if we prove that W4W3 is also trace class we 
can deduce from the cyclic property of determinants and (2.12) that 

(2.13) det(l-VKV^i) =det(l-W4W3). 

Recall from (2.9) that only the first column of (I + Wi)~iVP(T- + T+)(l + T+)-^K° has 
non-zero entries. Since U(l -|- T"'')PV~^ is upper triangular and is diagonal, the same 
is true for W4W3. Observe that V"^(l + Wi)-^V = {\ + PT+P)"\ so all the V's cancel in 
W4W3. For the first column of this operator-valued matirx we get using (2.9) that 

(W4W3)^^ = C/e"(*^-^*")'^P,„e-(*""*i)^5oC/^^ 

n-1 
i=k 

n-1 
i=k 



n-1 
i=k 



Telescoping the last sum yields 
(2.14) 

(W4W3)^ 1 = Ue-^'^-'^^^BoU-' - C/P,.,e-(*'=-*'=+i)^P,,^, • • • P,„e-^'"-''^^ BoU-' 



U 



-(ife-i„)A 



Using this last decomposition we get directly from the proof of Proposition 2.3(a) that 
(y\/4\N^) 1^ is trace class. This justifies the identity (2.13), and then since only the first 
column of W4W3 is non-zero we deduce that 

det(l - VKV-^) = det(/ - (W4W3)^ ^ 



L2(IR) 

The result now follows from the above formula for (\N4\N3j^ ^ with k = 1. □ 

Remark 2.2. A complete proof for the Airy2 case can be obtained from the above argu- 
ment by replacing — A by i?, by K\i, and both V and U by I. As we mentioned in 
Remark 2.1, this case presents the additional issue that the operators e*''^ involved in T"*" 
and T~ do not even map L^(M) to itself (in fact, note that H has the whole real line as 
its spectrum). T~, which is associated to operators e*^ with t > 0, presents no difficulty 
in the above proof. In fact, it always appears applied after K, which in this case is the 
diagonal matrix with K^i in each diagonal entry, so that since K\[ projects onto the nega- 
tive eigenspace of H (see Remark 1.1), each entry in T~K is a bounded operator acting on 
L^(M). This is analogous to the fact that, in the Airyi case, the operators e~*^ for t > 
always appear after Bq. 
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To deal with T"*" we start with the formula 

/OO fOO 
dy dAe^*Ai(x + A)Ai(y + A)/(y). 
-OO J — OO 

One can check that for any / G L^(M) the integral is convergent, and thus e~*^/ is well 
defined, though not necessarily in L^(M). The key is to notice, again using the formula, 
that for any a the operators Pae~^^ and e~^^ Pa are Hilbert- Schmidt (see (3.10)), so that 
Pae~^^ Pa = {Pae~^^){e~2^ p^'j is trace class by (2.18). In particular, this implies that 
the operator Wi defined in (2.3) (with V = I) is trace class in the Airy2 case. To make 
sense of (I + Wi)~^W2, as needed in (2.2), we can use (2.8) directly together with (2.7) to 
write 

n—i 

(2.16) [(|+Wi)-^W2]^_^. = ^^(-1)'= Yl ^x.e-(*'"*"i)^P,,^e-(*"i-*"2)^ 

k=0 i<ai<---<ak<n 

'-'■k—1 '^k 3 

(cf. (2.10)), where the same argument can be applied to show that each term is well defined 
and is in fact trace class. This allows to derive (2.12), and it is easy to check that deriving 
(2.13) via the cyclic property of determinants involves no new difficulties. 

A final remark is in order. The operator P^.^e^(*i"*2)^ • • • e~(*""i^*"-'^-PE^ appearing 
in (1.6) is ill-defined because, unlike in the preceding discussion, an operator of the form 
Pae~^^ Ph does not map L^(M) to itself. Hence (1.6) should be understood as a shorthand 
notation for 

HA2{h) <Xi,...,A2{tn) <Xn) 

(n n—i 
i=l k=0 i<ai<---<ak<n 



/ L2(M) 

which is obtained from the above proof by working directly with (2.8) instead of (2.9). 
Alternatively, one can rewrite 



lP(^2(tl) < Xl,...,A2{tn) < Xn) 



det I 



Jtl-tn)H _ p (tl-t2)Hp Jt2-t3)H_3 

L J /L'^( , 

The product inside this last determinant was shown to be trace class in Proposition 3.2 of 
[CQRll] (cf. Proposition 2.3 below). 

Going back to the Airyi process, we turn next to proving the existence of trace class 
operators which are conjugate to the ones appearing in (1.7) and (1.14). Given x = 
{xi, . . . , Xn) and t = (ti, . . . , t„) with tj < tj+i let 

(2.17) Af = P^,e^(*i^*^)^P^2e-(*2-t3)A . . . g-(i„.i-i„)Ap^^_ 

For the case ti = £ + ^^{r — £), i = I, . . . ,n, and Xi = g{ti) for some g G H^{[£,r]) we 
write 

Let U be the operator defined by Uf{x) = e~'^^'^~^^^ f{x). 

Proposition 2.3. Fix £ < r and let g G H^{[£,r]). 

(a) U{Bo - Afe-^^"-^^'^^Bo)U~^ and U [Bq - ^^e-^''-^^^BQ)U-^ are trace class op- 
erators on L^(R). 
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(b) \\U[Bo — ^^e~^^'~^'^^BQ^U~^\\^ is bounded uniformly in n. 

(c) Let Uk = 2^'. Then 

^lim \\U{Bo - Km'-^'-'^^Bo)U~' - U{Bo - A^e,rf-^^-'^^ Bo)U-\ = 0. 

Proof. The proof is similar to that of Proposition 3.2 of [CQRll], although here using the 
conjugated kernels is crucial. 

Assume first that g{t) = and write s = r — i. We begin by considering the second 
operator in (a). Let ip{z) = \J\^ and write 

Y{x,z) = (e^^-Af,,.j)(:E,z)e-2-(^(z)e-2- and W{z,y) = {e-'^Bo){z,yMz)-^e^''e^y\ 
Then 

u(^Bo - Af^ ,,]e-^^5o)t/^' = VW. 

Since 

(2.18) II^W^IIi < ll^lbllW^lb 

(see [Sim05]) it is enough to prove that ||F||2 < oo and ||W^||2 < oo. 
The estimate for ||W||2 is simple: using (1.9), 

„-4s3/3+2(x+y)s p -4:S^/3+2ys 

2_ / A:/^ , „. , „2n2 _ / , 2n2 



\W\\i= dxdy --^5 Ai(x + y + s^) = / dxdy — Ai{y + s 



/oo 
dye~'^''l^^^y^M{y + s'f. 
-oo 

The last integral is finite thanks to the bounds 

(2.19) |Ai(z)| < Ce-i^'^' for z > 0, |Ai(z)| < C for z < 

for some constant C > (see (10.4.59-60) in [AS64]), and thus ||Ty||2 < oo. 

For V, recalling that we are taking g{t) = 0, we may shift time by —{^ + r)/2 in the 
definition of A^^^j to deduce that A^^ ,^j = ^[_5/2 s/2]' ^"^^ then by (1.15) we have 

= /I- n(-./2)=.,K./2)=AK0 < on [-V2, V2]j . 



/47rs 

Therefore 

g-(x-s/)2/4s-2(x+j/)s , 

y(x,y) = ^{y) -^== n(-./2)=..,f,(./2)=,(Ki) > for some I G [-5/2, s/2] 

The last crossing probability equals e~^yl^ if x < 0, y < and 1 otherwise (see page 67 in 
[BS02]), and thus 



(2.20) \Vf^ = ^l (ixdy(l + y2)[e- 

+ / dxdy{\^y^){e- 



|2 ^ _ / ^-r- W„ n i r,-(x-j/)2/4s-2(x+j;)s] 2 

/m2\(-oo,0]2 

-(2:+J/)^/4s-2(x+3/)sl 2 

'(-oo,0]2 

Both Gaussian integrals can be easily seen to be finite, so we have shown that ||F||2 < oo. 

For the discrete time kernel we can use the same argument. To simplify notation we will 
write the proof for the kernel of the form A^ (with (7 = 0), the same proof works for 
A^. We decompose the kernel as 
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where 

Vn{x,y) = ifiy)- -^== n"(-./2)=x,S"(./2)=,(^"(«) < [-«/2, ^2]) 

and is a discrete time random walk with Gaussian jumps with mean and variance 
s/n, started at time —s/2 at x, conditioned to hit y at time s/2, and jumping at times 
= — s/2 + -^^s, i > 1 (in the case of a kernel this random walk is not time- 
homogeneous, but this does not introduce any issues below). We deduce that 



(e-(^-m _ A^' .){x,y) = ^e-(-y)V4s-2(x+,)s 




■ 6"(-s/2)=x,fe"(s/2)=3/ \ 

A simple coupling argument (see the next paragraph) shows that the last probability is less 
than the corresponding one for the Brownian bridge, and thus we obtain for ||e~('"~^)^ — 
A^ j,j||2 the same bound as the one we get for ||e-('^-^)^-A}'^^,j||2 from (2.20). This bound 
is, in particular, independent of n, so we have proved (a) and (b). 

To prove (c) we use again the above decompositions into VW and VnW. Our goal is to 
show that \\Vn^W - VW\\i ^ as /c ^ oo. Since \\Vni,W - VW\\i < HK^ - V^lbll^^lb by 
(2.18) and we already know that ||W^||2 < oo, all that is left is to show that 

UK, -^||2 ^0. 

k—^oo 

Couple the Brownian bridge b and the conditioned random walk by simply letting 
6"'=(t"'°) = h{t^^) for each i = l,...,nfc. Since the Brownian bridge hits the positive 
half- line whenever the conditioned random walk does, it is clear that 

^-(x-yf /As-2{x+y)s 

(2.21) \Vn,{x,y)-V{x,y)\ = = qn,{x,y), 



where qn^i^^y) is the probability that the Brownian bridge h{t) hits the positive half-line 
for t G [—s/2, s/2] but not for any t G {t^*", . . . , t^^}. Since every point is regular for 
one-dimensional Brownian motion, qn^{x,y) \ as A; — ?• oo for every fixed x,y, and thus 
the monotone convergence theorem yields (2.21). 

To extend the result to G ( [£, r] ) we note that everything in the above argument 
deals with properties of a Brownian motion 6(s) killed at the positive half-line. In the 
general case we will have by (1.15) a Brownian motion h{s) killed at the boundary g[s) or, 
equivalently, a process h{s) = h{s)—g{s) killed at the positive half-line. Using the Cameron- 
Martin-Girsanov Theorem we can rewrite the probabilities for h{s) in terms of probabilities 
for h{s). Since g{s) is a deterministic function in ff^([£, r]), the Radon-Nikodym derivative 
of h{s) with respect to h{s) has finite second moment, and thus by using the Cauchy- 
Schwarz inequality we get (a) and (b) from the above arguments. The convergence in (c) 
follows as well from the above arguments because it only depends on almost sure properties 
of the corresponding Brownian motion. □ 



3. Regularity and continuum statistics 



We now use the Kolmogorov continuity criterion to prove the Holder continuity of the 
Airyi process (we will explain later how to adapt the proof to the Airy2 case). An important 
technical problem is that the kernel appearing inside the determinant in (1.7) is not trace 
class. 

To apply the Kolmogorov criterion we have to get an appropriate bound on 
det(/ -Bo + Pae'^Pbe-'^Bo) - det(/ - Bq + PaBo). 
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To deal with the fact that the kernels above are not trace class, we have to conjugate by a 
kernel U as in Proposition 2.3. The resulting bound in terms of trace norms gets bad as 
a, 6 — >• — oo. To get around this, we use the Kolmogorov criterion in the following unusual 
form. 

Given a stochastic process X{t) and M > we denote by X^'^(t) the truncated process 

X'^\t) = X{t)l\x{t)\<M + ^'^X(t)>M - Mlx{t)<-M- 

Lemma 3.1. Let X{t) be a real valued stochastic process defined for t in some interval 
/CM. Assume that the following two conditions hold: 

1. There is a dense subset J of I such that \\mK^oo^{\X{t)\ <K \/t ^ J) = 1. 

2. There are a, /? > satisfying the following: for each M > there is an e > and 
c > such that 

E(|X^^(t)-X^(s)r) <c|t-s|i+^ 
for all s,t £ I with \t — s\ < e. 
Then X{t) has a version on I with Holder continuous paths with exponent ^. 

The lemma follows immediately from the usual Kolmogorov criterion, which, applied 
to 2, shows that there is a version of X{t) such that, for each M > 0, X^\t) is Holder 
continuous with exponent ^. Such a function cannot be discontinuous if it is bounded on 
a dense set. 

In view of this lemma, after we verify the first condition (which we do in the next result) 
it will be enough to consider the truncated process A^{t). Throughout this section all 
Fredholm determinants will be computed on L^(R), while c and d will denote positive 
constants whose values may change from line to line. 

Lemma 3.2. Fix L > and write Diin) = {^L, k = -2", . . . , 2"}. Then 

lim F{Ai{t) <M yte Un>oDLin)) = 1. 

M-5-00 

Proof. By Theorem 1, Proposition 2.3(c) and the bound 

(3.1) |det(/ + Qi) - det(/ + Qa)! < \\Qi - Qalliell'^illi+ilQ^Hi+i 
for trace class operators Qi and Q2 (see [Sim05]), we have 

F{Ai{t) <M yt£ Lln>oDL{n)) = lim F{Ai{t) <Myte Diin)) 

= det(/ -Bo + ^flL/2,L/2f~''''Bo), 

where Aj^^^^g l/2] denotes A^_^/2 L/2] '^i^h g{t) = M and, we recall, the operator inside the 
determinant is trace class after conjugating by U as in Proposition 2.3. Using (3.1) again 
we deduce that it is enough to show that 

(3.2) lim \\U{Bo-Al[L/2L/2f~''^Bo)U-\ = 0. 

Following the proof of Proposition 2.3(a) we have 

\\U{Bo - AflL/2,L/2f^''^Bo)U-\ < \\Vh\\Wh 

with V and W as in that proof. Recall that W does not depend on M and has finite 
Hilbert-Schmidt norm, so all we need is to show that ||T^||2 0. To estimate this last 
norm we can proceed exactly as in the arguments leading to (2.20), only replacing s by L 
and the barrier at for the Brownian bridge by a barrier at M, so that the corresponding 
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crossing probability is now e~^^~^^^^~^^'^^^ for x,y < M and 1 otherwise. We obtain, after 
some simple manipulations, 



+ ^ [ (ixdy(l + y2)re-(-+J/)V4L-2(x+y)L-2MLi2_ 
4ttL y(_oo,ol2 



'{-oo,0]2 

The last two integrals are easily seen to go to as M — )• oo, and (3.2) follows. □ 

Next we verify the second condition in Lemma 3.1. By the stationarity of Ai we may 
take s = 0. 

Lemma 3.3. Fix 6 > 0. Then there is a to G (0, 1) and no G N such that for < t < to, 
n > no and M = (31og(t~(^+")))^^'^ we have 

where the constant c> is independent of 6, no and to- 
Proof. By the stationarity of the Airyi process 



E([<(t)-<(0)]^"l^M(o)^^M(,)<_,,) <(2M)2"2P(A(0)<-M). 

Now P(^i(0) < -M) = Fgoe(-2M) < as M ^ oo by the results of [BBD08]. 

Hence we get 

E([<(t) - <(0)]^"1^m(o),^m(,)<_m) < c(2M)2"ti^-" < cti+(i-^)" 
if t is small enough. Thus it will be enough to prove the estimate 
(3.3) q{t) := E([<(t) - <(0)] '"1^M(o)^^m < ct'+^'-'^- 

for small enough t. 

Let F{a,b) = P(^i(0) < a,Ai{t) < b) and G{a) = P(^i(0) < b). Since ^G{aAb) = 
except when a = b we have 

/oo POO pp. 

da db{a-bf^——[F{a,b)-G{aAb)]. 
-M J-M dadb 

Truncating the upper limits at > for a moment and integrating by parts the integral 
becomes 

r da ({a - K)^^^[F{a,K) - G{a)] - (a + M)2"|-[F(a, -M) - Gi-M)]] 
J-M \ Oa oa J 

+ / da db2n{a-bf''~^ — [F{a,b)-G{aAb)] 
J-M J-M oa 

= -2 / da (2n(a - i^)^"-^ [F(a, K) - G{a)] - 2n{a + M)^"-! [F(a, -M) - G(-M)]) 

J-M 



K r-K 

da / db 2n(2n - l)(a - 6)2("-i) [F(a, b) - G{a A b)] 

J-M 
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(note that we have cancehed some boundary terms). We will see below in (3.9) that 

POO 

\F{a,K) - G{a)\ < cM^l^e'+'"'"" / dxe-^'/^ 

Jt-^/^{K~a) 

whence it is easy to see that the first integral on the right side above vanishes as K ^ oo. 
We deduce then that 

/oo 
da {a + Mf""-^ [G{-M) - F{a, -M)] 
-M 

/oo /"OO 
da / d6(a-6)2("-i)[G(a A6) -F(a,6)]. 

We will estimate the last double integral, the first integral in the last line can be estimated 
similarly. Since the integrand is symmetric, it will be enough to restrict the integral to the 
case — M < a < b. Using the definitions of F and G and Theorem 1 we have 

(3.5) F(a, b) - G{a A 6) = det(/ - Bq + PaC^^ PbC^*^ Bq) - det(/ - Bq + PaBo). 

Recall that the operator inside the first determinant is trace class after conjugating by the 
kernel U introduced in Proposition 2.3. We will use the bound 

(3.6) |det(/ + Qi)-det(/ + Q2)| < ||Qi - Q2||ie"'5^-'52"^+''l'^^"^+\ 

which follows directly from (3.1), to estimate the difference of determinants in (3.5), so our 
first task will be to estimate the trace norms of the operators 

Q2-Qi = U{Pae'^Pte-'^Bo-PaBo)U-^ and Qi = U{PaBo - Bop-^ 
for -M <a<b. 

We will use a different approach, and in particular a different choice of the kernel U, 
than the one used in the proof of Proposition 2.3. In what follows we will write x = 2^/^x 
and y = Let 

f//(x) = e-(*+-)^</.(x), where </'(x) = e^n^^.^vaM + 1.<-2V3m 
and a = M~^. We bound first the norm of Qi. Using the identity 

/oo 
du M[a + u)Ai{b - u) = 2~^/^ k\{2^^'^{a + b)) 
-oo 

we have 

= -2^I^Q\Ql with Q\{x, u) = l^>ae-(*+")*(^(x)-^ Ai(x + n)e(*+"/2)«, 

^^■^^ Ql{u, y) = e(*+°)j^<^(y) Ai(y - n)e-(*+"/2)-. 
Now (using the fact that a > — M) 

/"OO /"OO /"OO /"OO 

||Qi||2=/ dx due-2*^Ai(x + n)2e(2*+°)« = / dxe-^^*+°'^^ du M{uf e^'^^+''> . 

Ja J —oo Ja J—oo 

By (2.19) the last integral in u is bounded by c{t + a)""*^/^, and then 

HQ} lb < c{t + a)-3/4e-'=(*+")'^ < c'M^/\ 

where the second inequality follows from the choice a and M and the fact that a > —M. 
For Qi we have 



-oo 
oo 



\\Ql\\l= I dy dne2(*+")V(y)'Ai(y-u)2e-(2*+")'^ 

) J—oo 

POO 

dye'^^q^iyf / Ai(-u)2e-(2*+°)". 
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The u integral is bounded by c{t + a)~^/^ as before, while the y integral equals 

-M i-oo 



-M 

SO we also have ||(5i||2 ^ cAf^/^. Using these two estimates with (2.18) and (3.7) we 
conclude that 

(3.8) IIQi 111 < cM3/2. 

Now we need to bound ||C/(Q2 — Qi)t^ ""^lli- Recall that we are assuming a < 6, so that 
Paie'^Pb - Pbe'^) = -Pae'^Pb- Then 

/•OO -I 



UiQ2 - Qi)U-\x, y) = -l,<„e-(*+")^>(x)-i / dz ^ 

Jb Vint 

■ e-2*V3-(^+y)t Ai(^ + y + t2)e(i+a)S(^(y) 

/OO 1 
_ dz ^e-"/^l^,>,_,l.<.e-(*+")^<^(5)"i 

where we performed the change of variables z = x + \ftz. We regard this as an average of 
the kernels Cz{x,y) given by 

C,{x, y) = l^,>,_^^^<X£)-ie-2*'/3-(-+-)*-(^-S)*-"(--5)+*'^'^Ai(x + y + Viz + t^)m, 
SO that 

/OO 1 roc -I 

dz^e-'"/^C,\\i< I d~z^e-''"/^\\C,\\i, 
-OO V47r ^J^'K 

where the second inequality follows from the fact that Ci vanishes for \ftz < b — a. The 
same argument as the one used to estimate ||Qi||i with only a bit of extra arithmetic gives 
the same bound for ||Ci||i and thus we get 

\\U{Q2 - Qi)f/-ii < cAf3/2$(t-i/2(5 _ «)) 

with $(x) = dz e~^^ (in fact a better bound can be obtained in this case without 
much difficulty, but we will not need it below). 

Using the bounds on ||[/QiC/-i||^ and \\U{Q2 - Qi)U''^\\^ in (3.5) and (3.6) we deduce 
that 

(3.9) |F(a, h) - G(a A 6) I < cM'''^'^{t-^l^{b - a))e^^''^^'"'' < ct-^^r^/\b - a)) 
by our choice of M. Therefore 



OO ;>oo 



da / dbia- 6)2("-i) [G{a A 6) - F{a, b)] 

M J-M 

poo poo 

<ct-^ da d6(a-6)2("-i)$(t-i/2(5_a)) 



M J-M 
OO roo 



= ct"-M da dbia-bf^^'-^^^ib-a). 
J-M J-M 

Using the standard estimate <I>(x) < ce~^^/^ as x — )• oo it is not hard to see that the last 
integral is bounded by cM2("'~^). Using this in the second integral in (3.4), and recalling 
that a similar estimate holds for the first integral, we deduce that 

q{t) < cn2M2("-i)p-3 

and thus, using our choice of M, (3.3) follows. □ 
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Proof of Theorem 2. The last two lemmas allow to check the hypotheses of Lemma 3.1, 
which yields the result for the Airyi case. 

The proof for the Airy2 case is slightly simpler because the operators involved are trace 
class, and can be obtained by adapting the preceding arguments as we explain next. 

The one-point marginal of A2, which is given by the Tracy- Widom GUE distribution, 
satisfies the tail estimate FGUE(-Af) < ce~T2\^'\'^ (see [TW94]). Choosing now M = 

(l21og(t~''"'^^"^))"'^^"^ it is not hard to check that the main argument used in the case of the 
Airyi process works in exactly the same way if we change our determinantal formulas to 
the corresponding ones for A2- Thus all we need to do is to obtain an analogous estimate 
on the difference 

F{a, b) - G{a Ab) = det(/ - Km + Pae-'^ Pbc'^ Km) - det(/ - Km + Pai^Ai) 

for —M < a <b. Recall that the operators inside these determinants are trace class in this 
case, so there will be no need to conjugate. Proceeding as in the proof for ^1 we need to 
bound the trace norms of the operators 

Q2-Qi = Pae-'^ Pbc'^ Km - PaKM and Qi = PaKM - Km- 

We start with Qi, which we rewrite as —{Pae^^^^ N){N^^e°'^ Kmi) with a = and 
N the multiplication operator Nf(x) = ip{x)f{x) with ip{x) = (1 -|- x^)^/^ (the choice of ip 
is not particularly important). It is easy to check (see (3.3) in [CQRll]) that 

||iV-^e"^KAi||2 < ca"^ 
for some c > 0. On the other hand, 
(3.10) 

poo roo p 

||p„e-"^||2 = / dx dy dA dA e^"(^+^) Ai(x - A) Ai(y - A) Ai(x - A) Ai(y - A) 

poo rco roo poo 

= dx (iAe-2"^Ai(x-A)2= / dxe-^^M dAe'^"^ Ai(-A)2 

Ja J —00 Ja J —00 

-3/2 -2aa 

where we used (2.19) as before. Using these two bounds together with (2.18), our choice 
of a and the fact that a > —M, we get 

(3.11) IIQilli < ca-^/^e"'^'^ < c'M^/\ 

We turn next to the trace norm of Q2 — Qi- Recalling that H = —A + x and defining 
the multiplication operator (e"^/)(x) = e"^/(x) (the reason we use the letter ^ instead of 
X in the definition is that we will use the operator at points other than x below), one can 
derive formally, using the Baker-Campbell-Hausdorff formula, that 

where e^^^ f{x) = f{x + t^) (see [QRllb] for a similar computation). This formula can 
then be checked directly by integration using (2.15) and therefore we may write, similarly 
to the Airyi case, 

/OO -1 
-00 V47rt 

/OO 1 
d~z -^e-^"/^ (e*^/3+*^^e-*«ne*^KAi) {^ft~z + x, y) 
-00 v4vr 

1 -2 
dz—=e~'' ''^Cz(x,y), 
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where Cj = Pae*^/'^"'"'^^^"^*^)^e^*^Pf,e*'^ilAi and we have used the fact that Cz vanishes for 
^/tz < h — a — t^. Proceeding as above we write, with a = M~^, 

where || • ||op denotes the operator norm in L^(M) and we have used (2.11). The first norm 
on the second line can be easily bounded by ce~^*^/^~*''~*^''^^, while for the other two norms 
we have already obtained ||Pfee(*-°)^||2 < c(a - t)-3/4e(*-°)'' and ||e"^KAi||2 = {2a)-^/^ 
in the derivation of (3.11). Since we are only interested in the case ^/iz > b — a — t'^, we 
have e-^tys-tb-t^/^z < ^ty3-2tb+ta ^^^^ 

lla-lll < C{a - t)-3/4^-l/2//3-(t+a)fe+ta < ^/^5/4^ 

where we have used the again our choice of M and a and the fact that — M < a < b. 
Plugging this in the above formula for Q2 — Qi we get 

IIQ2 - Qilli < cA/5/^$(t-^/2(6 -a- t^)). 

This estimate, together with the one for ||(5i||i, allows to derive the an estimate analogous 
to (3.9): 

\F{a, b) -G{a^b)\< cM^/^^{r^/^{b -a- t2))ei+^*^'^' < cr^<^{t-^''^{b -a- t^)). 

Comparing with (3.9), the only difference is the additional shift by —t^/"^ in the error 
function <1>, but it is easy to see that this does not introduce any difficulty, and the rest of 
the proof follows as for ^1. □ 

Finally we turn to the continuum statistics formula for the Airyi process. 

Proof of Theorem 4- With the notation introduced before Proposition 2.3 we have 

nMh) < gih), . . .,Ai{tn,) < 9{tn,)) = det(/ - Bo + A^^^[,^,]e-(^-^)^Bo)^^^^^ , 

where = 2^. Since, by Theorem 2, Ai has a continuous version, the probability on the 
left side converges to F{Ai{t) < g{t) \/t € t]), and thus it is enough to show that 

limdet(/-[/{B„-A»,,„,.-<'-')-B„)C/-)^^^^^ 

= det(/-C/(B„-Af,,,.-<-')-B„)[/-')^,^_,^, 

where = 2^ . Since A 1— ?• det(/ + A) is a continuous function on the space of trace class 
operators by (3.1), the identity follows readily from Proposition 2.3(c). □ 

4. Local Brownian property of Airyi 

Note that, by stationarity and time reversibility, it is enough to study the finite dimen- 
sional distribution of Ai at times s = < ii < • • • < t„. We have the following formula for 
the Airyi process conditioned at a point. 

Lemma 4.1. For < ti < • • • < i„, 

(4.1) ¥{Ai{ti) < X + yi, . . . , Ai{tn) < x + y„ | ^(0) = x) 

^ P(^i(0) < x,Ai{ti) <x + yi,...,Ai{tn) <x + yn) 



^-^goe(^^) 



tr 



1 



I-B^ + A|jf+")e-*"^i?o) 5,e*^^At^+"e-*"^i?o 
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where is defined in (2.17) and (0,t) and (x,y + x) are notations for the vectors 
{0,ti, ...,tn) and {x,yi+x,...,yn + x). 

Note again that the analogous formula is true for Airy2. We remark that in the trace ap- 
pearing in (4.1) we should be conjugating by the operator U introduced before Proposition 
2.3 to make sure that the operator is trace class. The same is true for the calculations that 
follow. To simplify the argument we will ignore these conjugations and skip some details 
throughout this section, we hope that at this point the reader can fill in the necessary 
arguments. 

Proof of Lemma 4-i- Note first that 

P(A(ti) <x + yi,...,Ai{tn) <x + yn\Ai{0) = x) 

GOE V / 



h=x 



where we have used the fact that P(^i(0) < x) = Fqoe{2x) and Theorem 1. Now recall (see 
[Sim05]) that if {^(/3)}^>o is family of trace class operators which is Frechet differentiable 
(in trace class norm) at (3 = h then 

(4.2) dh det(/ + A{h)) = det(/ + A{h)) tr[(/ + A{h))-^ dhA{h)]. 

The result now follows from computing the Frechet derivative of ^(^q\) \ which can be 
shown without difficulty (after introducing the necessary conjugations) to make sense in 
trace class norm. □ 



Proof of Theorem 3. We study the last line of (4.1) and to make it easier to read we call 



L = Bq + A[Q'^+^^e-*"^So- Note first of all that it is given exphcitly by 



tr 



(/ - L)-^ 5,e*i^A^+"e-*"^So" 

oo i-oo 

dz e*'^P^+y, ■ ■ ■ e(*"-*i)^P,+j,Jx, z) / due-'-'^Boiz, u) {I - L)-\u, x). 



-oo ^ — oo 

Shifting z hy X and using the translation invariance of the heat operators we can rewrite 
the trace as 

/oo fOO 
dze^i^Pj^i •••e(*"-*"-i)^Pj^„(0,z) / due-^"^Bo{z + x,u) {I - L)-^{u,x). 
-oo J — oo 

If we put in the Brownian scaling t i— )■ et, y i— )• i/ey we get 

/oo roo 
dz e^*^^P^,^ • • • e^(*"-*i)^P^j^„ (0, z) / du e~''"'^Bo{z + x, u) {I - L,)-\u, x), 
-oo J — oo 

where is defined in the obvious way by introducing the Brownian scaling in L. Since 
the heat operators are invariant under this scaling we can change z i— )■ y/ez to see that this 
is equal to 



/oo POC 
dze^i^Pj,, •••e(*"-*"-i)^P3,„(0,z) / due-''"'^Bo{^z + x,u) {I - Ley\i 
-oo J — oo 



fOO 

, , . (u,x]. 

-oo 

Combined with ^Fgoe(2x) = -Fgoe(22;) du Bq{x, u) {I - Bq + PxBq) ^(n, x), which 
follows easily from (4.2), we obtain (1.10) from this and (4.1). Now (1.12) goes to 1 as 
e — )• by the continuity of Airyi proved in Theorem 1. On the other hand, one can show 
that converges to L as e — t- in trace class norm, which implies (see [Sim05]) that 
(/ — Lir)~^ —)•(/ — L)~^ in the same sense. Using this it is not hard to show by the domi- 
nated convergence theorem that (1-11) goes to 1 as e — )• 0. This implies the convergence of 
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the finite dimensional distributions to those of Brownian motion, and thus concludes the 
proof. □ 
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